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Abstract 

We extend form-factor perturbation theory to non-integrable deformations of mass¬ 
less integrable models, in order to address the problem of mass generation in such 
systems. With respect to the standard renormalisation group analysis this approach 
is more suitable for studying the particle content of the perturbed theory. Analo¬ 
gously to the massive case, interesting information can be obtained already at first 
order, such as the identification of the operators which create a mass gap and those 
which induce the confinement of the massless particles in the perturbed theory. 


Introduction. Given the large number of remarkable results obtained from the study 
of two-dimensional integrable quantum held theories (IQFTs), at present one of the most 
interesting challenges consists of developing a systematic approach to study non-integrable 
models, at least when they are deformations of integrable ones. For massive held theories 
a convenient perturbative scheme, based on the exact knowledge of the form-factors (FFs) 
of the original integrable theory, was suggested in P] . Already at hrst order, it proved able 
to provide a great deal of information, such as the evolution of the particle content, the 
variation of the masses and the change of the ground state energy - results successfully 
checked by numerical studies. 

The main purpose of this paper is to extend Form Factor Perturbation Theory (FFPT) 
to non-integrable deformations of massless IQFTs. The most fundamental question that 
one may ask in this context is whether a perturbation creates a gap in the excitation 
spectrum - a problem usually addressed via the renormalisation group (RG) equations 
near a hxed point [2]. Moreover, if massive particles are created, one would like to 
understand whether they are adiabatically related to the original massless excitations or, 
like in the massive case, confinement takes place. Since the RG eq.s cannot provide a 
complete answer to any of the above questions, it is worth exploring other alternative 
routes. The FFPT relies directly on the particle description of the unperturbed theory 
and, for this reason, it seems to be the most natural and suitable approach for studying 
the evolution of the particle content when the perturbation is switched on. 

Our analysis is presently limited to the first order of FFPT, its extension to higher 
orders being, as in the massive case, an interesting but non-trivial mathematical prob¬ 
lem. Despite the fact that one must be careful in handling results at such low order, 
some useful conclusions can nevertheless be reached. For instance, it will be possible to 
discriminate between operators which do not spoil the massless nature of the theory and 
those which instead induce a mass gap in the spectrum. Moreover, the conhnement of 
the original massless excitations can be traced back to the non-local properties of the 
perturbing operator with respect to them. These results provide the first information on 
the perturbed theory and may guide a further analysis of its properties. It should be 
stressed that answering the above questions in their full generality is, obviously, a fairly 
complicated problem since it concerns the global structure of the RG flows rather than 
their local properties around the fixed points. It is well known, for instance, that adding 
a relevant perturbation to a massless action does not necessarily imply that the resulting 
infrared theory will be massive: indeed the perturbing operator may induce a flow into a 
new critical point, with some of the massive excitations decoupled from the new massless 
ones pin. An example even more subtle is given by the roaming trajectories discovered 
by Ah Zamolodchikov [H] (and further analysed in pm), i.e. an inhnite cascade of 
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massless flows finally ending in a massive phase. 

FFPT for Massive Field Theories. Consider non-integrable theories obtained as 
a deformation of an integrable action Aint 

A = Ant + A j (fx'^ix) . (1) 

The exact knowledge of the FFs of the operator \l/(a:) on the asymptotic states of the 
unperturbed theory allows one to set up an expansion of various physical quantities of 
the new theory in powers of A, the so called FFPT p. We present initially some known 
results of FFPT for massive theories in a way that is more suitable for the extension 
to the massless case. It is useful to recall that in most of the cases of interest, the 
integrable action of a massive theory can be defined in terms of a deformation of a CFT 
|Hj, Aint = AcFT + g f (Px <h(x), where <h(x) is a relevant scalar held of conformal weights 
A$ = < 1. 

Let us hrst assume that the theory has only one massive particle in the spectrum, d(/3), 
where [3 parameterises the dispersion relation: pq = m cosh (3 and pi = m sinh j3. The 
integrability of the theory allows one to compute its exact factorised scattering amplitudes 
0, *S'(/9 i2) (/9i2 = Pi — P 2 ), and the FFs P of the various operators O on the set of 

asymptotic states 

F^(A, /?2,..., Pn) = (0 I O{0) I A{P,)AiP2) • • • A{Pn)) . (2) 

A convenient way to study the mass correction induced by the non-integrable deforma¬ 
tion CD is to employ the Hamiltonian formalism, in the same spirit of standard quantum 
mechanics perturbation theory. The Hamiltonian associated to 0 can be written as 

H = -^ J Too(x^, 0) 0) , (3) 

where T^y{x) is the stress-energy tensor of the integrable theory, Aint- The operator Tqo 
can be expressed in terms of the its trace, ©(x), using the conservation law = 0, 

- do) Too(x\ a;°) . (4) 

In particular, for the two-particle Form Factor we have 

(0 I Too(x\a:°) | A{pi)A{pj)) = - sinh^ ^ (0 | 0(a:\x°) | A{pi)A{pj)) . (5) 

The essential results of FFPT are easily re-derived within this formalism. Let us first 
consider the unperturbed integrable case, A = 0. Evaluating the matrix element of both 
sides of Eq. (jS]) on the asymptotic states {A{Pi) \ and | A[Pj)), and using the relation ©, 
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one obtains the usual normalisation condition for the trace of the stress-energy tensor of 
the massive integrable theory, 

{A{f3) I 0 ( 0 ) I A{j3)) = F®(z7r) = 27rm^ , (6) 

an equation which shows the relationship between the FF of this operator and the mass 
scales of the theory. 

Repeating the same procedure for the non-integrable theory ©, one obtains instead 
the first order correction to the mass of the particles, as given in pQ, 

~ 2XF^{iTi) . ( 7 ) 

If the operator \h(a:) is non-local with respect to the particles ^(/9), F^{I3) has a pole for 
(3 = iiT and © diverges. This divergent correction to their masses implies the conhnement 
of the particles A(/9), that are no longer excitations of the action ([Q) p. This phenomenon 
appears for instance in the magnetic deformation of the low-temperature phase of the Ising 
model nmniini as well as in two-frequency sine-Gordon model jl2j . 

It should be noticed that, since this is a strong coupling analysis, i.e. carried out in 
the infrared region (IR), if A in ([1} scales under RG, it has to be replaced in 0 by its 
renormalised value at energy of the order of the mass of the theory 

A ^ A"^^ ~ A(m-i) . (8) 

As a consequence, unless the RG flow is known exactly, quantitative predictions can be 
made only on universal mass ratios. 

If the theory has n non-degenerate particles, Aa(/3), with masses nia (a = 1,... ,n; 
'ma 7^ 'ma') the above analysis can be easily extended and gives the following mass variation 

8ml ~ 2XFl{i7r) . (9) 

When some of the particles, say n' {p! < n ), have the same mass, this equation has to be 
generalised like in quantum mechanics perturbation theory for degenerate levels, i.e. the 
perturbed masses are obtained by diagonalising the matrix where 

indices k, I belong to the degenerate multiplet. If the symmetry of the perturbing operator 
is less than the symmetry of the multiplet, the perturbation will typically split it. 

Massless IQFTs. Massless non-scale invariant IQFTs are associated to RG flows 
between two different fixed points. With respect to their ultraviolet fixed point, such 
theories admit a well-defined description in terms of the corresponding GFT perturbed 
by a relevant operator. However, from the physical point of view of selecting the low- 
energy massless excitations, it is more appropriate to view them as irrelevant perturbation 
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of their IR fixed point action 


*^massless *^CFT 9 I d X -f- * • • 


/ 
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where the irrelevant held $(a;) specihes the approaching direction to the CFT of the IR 
hxed point. The scattering theory of massless IQFTs is discussed in detail in |1] whereas 


their Form Factors in PI, so we shall outline only some basic facts below. The excitations 
of these theories consist of right {R) and left (L) moving particles. They are dehned as 


Pi > 0 and Pi < 0 branches of the relativistic dispersion relation po = |pi |, which can be 


parameterised as Po = Pi = {M/2) e^ for the R movers, Afi{j3), and po = —Pi = {M/2) e~^ 
for the L movers, Al{/3), where M is a mass scale. Within this parameterisation, the 
Maldenstam variable for the RL scattering process is given by: SRL{l3ij) = M^e^*-’'. 
Contrary to the massive case, where the threshold of the scattering process is given by 
[3ij = 0, for the RL sector of the massless scattering the threshold is reached in the limit 
jSij —> —cx). In the RR and LL sectors the Mandelstam variable is always zero, showing 
that all analyticity arguments of the S'-matrix theory cannot be applied: the scattering 
amplitudes in these channels can be properly dehned only as analytic continuation of the 
massive case |1]. For this purpose, in fact, it is useful to regard the massless excitations 
as a particular limit of the massive particles^. 

Writing the S'-matrix in a compact form 74 ai(/ 5 i)Aa 2 (/ 52 ) = Sa^^a 2 {Pi 2 ) ^Q! 2 (/^ 2 )^ai(/dl); 

(ctj = R, L) the equations satished by the FFs can be written in analogy to the massive 
case Pl. For the two-particle matrix element {/3 12 ) — (0|C>(0)|Aa^(/3i), Ra2(/32)), 

we have for instance 



( 11 ) 


where 70 is the non-locality index of the operator O with respect to the massless particles. 
Their analytic structure, however, differs from the massive case. In massive theories, the 
multi-particle FFs are meromorphic functions in the strip 0 < \m.(3 < 2 t[ and present 
simple pole singularities associated either to bound states or to particle-antiparticle an¬ 
nihilation processes. In massless theories the same kinds of singularities are expected in 
the RR and LL sectors, since they formally behave like the massive cases. In the RL and 

^For instance, if A(/3) is a massive excitation of mass m with a Wmatrix equal to S{(3), the massless 
limit is constructed by shifting the rapidities j3 —> Pr^l i /3o/2 and taking the double limits Po —>■ 00 
and m —> 0 while M = me^° is kept fixed: Ar^r{P) = \m\f}^^aoA{P ± /3o/2). When one considers the 
S-matrix in the RR and LL sectors, the rapidity shifts cancel and therefore Srr{P) = Srl{P) = S{P). 
As functions of the rapidity variable, these amplitudes are then expected to satisfy the same equations 
valid for the massive case. 
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LR sectors instead, bound state poles are absent while kinematic poles may appear only 
if both particles have vanishing momentum. However, instead of producing a recursive 
equation like in the massive case, here the presence of kinematic poles imposes a condition 
on the asymptotic behaviour of the FF’s. In particular, like in the massive case, a pole is 
present in the two-particle FF only if the operator is non-local 

lim + /3rl) = oo for non — local (12) 

PRL^-OO 

(an analogous equation can be written for the FF in the LR sector.) 

FFPT for Massless Field Theories. Suppose now that an operator ^(x) of the 
infrared CFT is added to the effective action m, so that its integrability is broken 

A= Amassless + >^ j d‘^X^{x) . (13) 

Repeating initially the analysis of the previous section for the unperturbed case A = 0, 
one hnds^ the following normalisation conditions for the trace of the stress-energy tensor 

= FUi^) = 2 ^ ; ( 14 ) 

FidH) = 0 . ( 15 ) 

The last equation can be viewed as an essential property of a massless integrable theory, 
i.e. a non-trivial generalisation to massless non-scale invariant theories of properties of 
CFTs. 

Consider now the case when A is non-zero. If the perturbing operator fli(a:) has van¬ 
ishing FF on the RL (LR) sector, it is easy to see that, at the lowest order, it does not 
change the masslessness nature of the theory. Indeed, it does not spoil both the validity of 
eq. m and the analytic structure of the Green’s functions of the original massless theory. 
On the contrary, if the operator (^(x) has non-vanishing FFs in the RL (LR) sector of 
the theory, this perturbation immediately generates a mass gap, a quantity which can be 
estimated by sandwiching Eq. (0 on R and L asymptotic states 

(5m ~ 2 lim ■ (16) 

Prl^-oo 

From a kinematical point of view, the above limit is the expectation value of the perturbing 
operator at the (zero-energy) threshold of the crossed RL channel. In the above equation 
the effective coupling constant, A®-^-^, is dehned like in (|H} with m —> 0. As a consequence, 

^It should be kept in mind that to avoid trivial vanishing of the RR (LL) FFs in taking the massless 
limit of m, one has to rescale scalar operators by their mass dimension 0{x) —> 0{x)/nF‘^° and define 
their FFs Fi(q,(/3i 2) = limm^o , where F^{j3i2) is the two-particle form-factor of the massive 

version of the theory. 
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if the perturbing operator is irrelevant with respect to the IR CFT, scales to zero, i.e. 
the actual mass gap vanishes in the infrared region although it is present at intermediate 
scales. On the other hand, if it is relevant, it will grow to the scale of the mass being 
generated. At this level the relevance of an operator has to be established by scaling 
arguments. In summary, two conditions have to be fulfilled for generating, at lowest 
order, a mass gap in the theory: the perturbing operator has to be relevant with non¬ 
vanishing RL (LR) matrix elements. 

Moreover, like in the massive case, the mass correction 5m may be a hnite or a di¬ 
vergent quantity, depending on the locality properties of the perturbing operator with 
respect to the fields that generate the massless particles. If the operator is local, 
5m is finite and the massive excitations of the perturbed theory are adiabatically related 
to the massless particles of the original one. If, instead, \l/(x) is a non-local operator, 
it follows from Eq. m and (US) that 5m diverges: in this case, the original massless 
excitations are confined as soon as A is switched on. In other words, the massive particles 
of the perturbed theory are not associated, in this case, to the operators that create the 
original massless ones. The examples discussed below should help in clarifying these two 
situations. 

If the theory contains more than one type of massless particle Aa,R/L (a = 1,... ,n) 
the previous approach has to be generalised in analogy with perturbation theory for de¬ 
generate levels. Since the particles Aa^a form a complete basis for the scattering theory, 
by using FFPT it should be possible, in principle, to predict whether any massless excita¬ 
tions survive, and this is a clear advantage with respect to the RG. Although a complete 
answer to this question involves the entire series in A, nevertheless the first order of FFPT 
may provide useful hints on the decoupling of massive and massless modes of the theory 
under investigation. 

Massless flows between minimal models. Let us now apply the above methods 
to some specific examples, starting from the massless flow between the Tricritical Ising 
Model (TIM) and the Critical Ising Model (CIM). The quantum field theory associated 
to this RG flow can either be seen as TIM perturbed by its sub-leading energy operator e' 
of conformal dimensions A^/ = A,./ = 3/5 or as CIM perturbed by the irrelevant operator 
TT (see j3] and references therein). The factorised scattering theory for this massless 
flow was first proposed in jl] and the basic FFs calculated in The spectrum consists 
of massless neutral fermions, with S-matrix SRPtiP) = Sll{(5) = —1, while Srl{(5) = 
tanh {(512 — iiT/A). As well as studying the non-integrable theory obtained by the insertion 
of the energy operator e(x) of the CIM, we will also consider the deformation of the 
massless action by the disorder operator ju(x). The latter is non-local, 7 ^ = 1/2, with 
respect to the massless fermion excitations. The energy operator has two particle FFs 
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only in the RL sector of the form na 


niiP) 



dt 

t 


. 9 

Sin 


((iga) 


sinh t cosh 



(17) 


where is a normalisation constant. The disorder operator, on the other side, has also 
FFs in the RR and LL sectors, of the form tanh(/5/2). For the FFs in the 

RL sector we find 


FrM = 



dt 

t 


sin^ 




sinh t cosh 



( 18 ) 


The above results agree with the roaming limit of the FFs of the sinh-Gordon theory na 
- a limit in which the sinh-Gordon model corresponds to the above massless flow laiisiizi. 

Using now (HED and it is easy to see that the perturbation by e{x) induces a 

hnite mass in the system, as it could have been expected on different grounds. At the 
critical point, in fact, e{x) is bilinear in the fermionic operators that generate the massless 
particles, e ~ 'ipip, and therefore the perturbed theory describes massive fermions in the 
presence of an irrelevant perturbation TT. 

Gonsider now the perturbation of the massless action by the non-local operator /i(x). 
By computing the limit (uni) of the two-particle FF (Hi of this operator, one sees that in 
this case 5m diverges, i.e. the initial excitations can no longer propagate as asymptotic 
states in the new vacuum of the theory created by the insertion of this field. Like in 
the massive case, there is a simple explanation of this conhnement phenomenon in terms 
of the LG effective description of the theory [T]. Indeed, in the unperturbed theory 
the elementary excitations can be equivalently considered as massless kinks interpolating 
between two degenerate minima of the LG potential ^3]. However the insertion of the 
disorder magnetic operator ^{x) lifts the degeneracy between the minima, thus making 
the kinks unstable. 

As a matter of fact, the flow between the TIM and the GIM is the simplest example of 
a one-parameter family of RG trajectories interpolating between the conformal minimal 
models Ap, with central charge Cp = 1 — G/p{p +1) (p = 3,4 describe the CIM and TIM 
respectively). The flows start from Aoo and pass close to all the other minimal models, 
remaining massless all the way down to the very last hxed point, p = 3, after which 
they become massive m ini nj. The trajectories going out from each critical point are 
described as Ap perturbed by the operator 0 i 3 jS]: Ap^^ = Ap + X J d^x^j’g (where the 
upper index in 013 indicates the relative GFT) and the excitations are massless kinks 
interpolating between the (p — 2 ) degenerate vacua of the effective LG potential ^3]- An 
interesting problem consists of predicting the evolution of the spectrum along these flows. 
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in particular the successive decoupling of the massless modes in the cascade of massless 
RG flows, by applying FFPT. The analysis of this problem is however beyond the scope 
of the present letter. 

Spinon confinement in sigma models. Another important application of FFPT is 
in the study of the mass spectrum of the 0(3) non-linear sigma-model with a topological 
term 


•Af) — 


2/2 


d‘^x + i9T, (a = 1, 2, 3 ; = 1) 


(19) 


where / and 6 are dimensionless coupling constants and T is the integer-valued topological 
term related to the instanton solutions of the model. The two values 6 = (0, tt) are the 
only ones for which the action (unD is known to be integrable. At 0 = 0 the excitations 
form a massive 0(3) triplet whose scattering theory was constructed in jSlE]. At 6 = tt 
the theory is instead massless [m ca eg 1201 and corresponds to the RG flow between 
the c = 2 GFT and the SU{2)i Wess-Zumino-Witten (WZW) model. The factorised 
scattering theory was suggested in irzi: it consists of right and left doublets, Aa^^i and 
Aa,L (a = 1, 2), that transform according to the s = 1/2 representation of SU (2) (spinons). 
However, as soon as one moves away from 6 = tt, the spinons conhne [2111221 and the 
actual spectrum of the theory in the vicinity of this point has been determined in [221 • 
Let us discuss in some detail how the spinon conhnement takes place. In the FFPT 
this amounts to show that the topological term is non-local wrt the helds that create 
the spinons, a property that can be easily checked by looking at the GFT limit of these 
operators. 

Gonsider A-jr as our unperturbed IQFT. Glose to the IR hxed point the massless 
flow can be described as a SU{2)i WZW model perturbed by the marginally irrelevant 
perturbation (Trgf)^ |T21 UHl CHI, A^^-^ = Asu{2)i + 1J d‘^x{TT (7 > 0), where g is 
the SU(2) matrix held. In terms of this formulation the perturbation that moves the 
topological term away from 6 = n is proportional to Tig ^Hl, i.e. to the only relevant 
SU (2) invariant operator in the theory that breaks parity. Thus, in the vicinity oi 6 = tt, 
the model is described by the effective action 


j^eff ^ ^ ^ 


d^xTr g , 


( 20 ) 


where 7 is a function^ of {6 — tt) that vanishes when 6 = tt. 

As discussed in Ref. the spinons are created by the primary operator 0^(^) 
(0^(^)) with scaling dimension h = (1/4,0) (h = (0,1/4)). They enter the operator 

^The form of this function determines the dependence on {9 — tt) of the mass gap m, since it scales 
as m ~ up to logarithmic correction m In a recent paper m it has been suggested that the gap 
behaves like {9 — Tr)^/^^ which would imply a dependence of g on {9 — tt) that it is not linear, as usually 
assumed. 



product expansion (OPE) 

0"(z)0^(ta) = {—y{z — (l + l/2(z — w)‘^T{w )...) 

-(-y(z - (r(w) + i/2(z - w)ar(w) ...), ( 21 ) 

where = —e~~^ = 1 , are the generators of the algebra, and q takes the values 

g = 0 for states that are created by an even number of spinons and g = 1 if the number 
of spinons is odd. The OPE between the spinon operator and the SU{2) currents J°‘{z) 
is standard: J°'{z)(j)'^{w) = {w)/{z — w) + ... . From these OPEs it follows that 

and 0“ are mutually local while 0“ and 0^ are not. In fact taking 0“(z) around 0^(tc) 
by sending z —>■ ze^” it produces a factor with 7 ^ = 1 / 2 . 

Therefore (Trg)^ ~ J“J“ is local with respect to the spinons and this explains why 
they are the fundamental excitations of A-„-, regardless of whether the perturbation is 
marginally relevant or irrelevant. Since Trg is proportional to (0’''0“ + 0“0^), the OPE 
m implies that this operator is instead non-local with respect to the spinons. Hence they 
get conhned as soon as the operator Trg is added to At,, i.e. the perturbed model (EOD 
has no longer spin 1/2 excitations. As discussed in j22], the actual massive excitations 
of the 0(3) sigma model with 6 *-term consists of a triplet of particles and a singlet, the 
former stable for all value of 9 whereas the latter stable only in an interval of values of 9 
near 9 = w. 

We are grateful to G. Delhno for important discussions and A. Nichols for reading 
the manuscript. D.C would also like to thank A. Tsvelik and K. Schoutens for useful 
discussions. This work is within the activity of the European Commission TMR program 
HPRN-CT-2002-00325 (EUCLID). 
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